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Abstract 

We derive the effective action for pseudoscalar mesons by integrating out vector and axial- 
vector collective fields in the generating functional of the bosonized NJL-model. The corre- 
sponding modifications of the nonlinear effective Lagrangian and the bosonized currents, arising 
at 0{p^), are discussed. 







1. Introduction 



A renewal of interest in chiral Lagrangian theory was excited by recent progress in the construction 
of reahstic effective chiral meson Lagrangians including higher order derivative terms as well as the 
gauge Wess-Zumino term from low-energy approximations of QCD. The program of bosonization of 
QCD, which was started about 20 years ago, in the strong sense is of course also beyond our present 
possibilities. Nevertheless there is some success related to the application of functional methods 
to approximate forms of QCD (see 0-0 and references therein) or to QCD-motivated effective 
quark models [0-|17 which are extensions of the well-known Nambu-Jona-Lasinio (NJL) model 
fl^ . These functional methods can be applied also to the bosonization of the effective four-quark 
nonleptonic weak and electromagnetic-weak interactions with strangeness change \A.S\ = 1 by using 
the generating fmictional for Green functions of quark currents introduced in jl^ , pof . 

The NJL model, which we consider in this paper, incorporates not only all relevant symmetries of 
the quark flavour dynamics of low-energy QCD, but also offers a simple scheme of the spontaneous 
breakdown of chiral symmetry arising from the explicit symmetry breaking terms due to the quark 
masses. In this scheme the current quarks transit into constituent ones due to the appearance 
of a nonvanishing quark condensate, and light composite pseudoscalar Nambu-Goldstone bosons 
emerge accompanied also by heavier dynamical vector and axial-vector mesons with correct relative 
weights arising from renormalization. The composite vector and axial-vector resonances are naturally 
required if one wants to describe the low energy aspects of QCD in a wider energy range up to typical 
masses of O(lGeV). (For other approaches to the problem of introducing chiral couplings for vector 
and axial-vector mesons to Goldstone bosons into the effective chiral Lagrangian, see for example 
refs.Q-Q and references therein.) 

Independently from the method of including the vector and axial-vector fields in the effective 
chiral Lagrangian, integrating out the heavy meson resonances essentially modifies the coupling con- 
stants of the pseudoscalar low-energy interactions. In particular, in rcfs.||2^, p3| it was shown that 
the structure constants Li of the Gasser-Leutwyler general expression for the 0{p*) pseudoscalar 
Lagrangian are largely saturated by the resonance exchange contributions giving a product of 
terms of 0{p^). But in this case, if the O(p^) Lagrangian contains meson resonances, their elimina- 
tion can lead to the double counting mentioned in ref . [^2| . The resonance contributions to the purely 
pseudoscalar chiral weak Lagrangian and the modification of its structure, induced by integrating 
out the heavy meson exchanges, were discussed recently in refs.|25 in both the frame of the 
factorization approximation and in the weak deformation approach. 

In this paper we consider the effective nonlinear Lagrangian for pseudoscalar mesons which arises 
after integrating out the explicit vector and axial-vector resonances in the generating functional of 
the bosonized NJL model. To perform such integration we use a method based on the invariance of 
the modulus of the quark determinant under a chiral transformations and on the application of the 
static equations of motion to a special configuration of the chiral rotated fields. The elimination of 
vector and axial-vector degrees of freedom from the modulus of the quark determinant leads to a 
modification of the general structure of the effective strong Lagrangian for the pseudoscalar sector 
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at O(p^) and to a redefinition of the corresponding Gasser-Leutwyler structure coefficients Li. This 
method of reduction of meson resonances can be extended to the procedure of chiral bosonization 
of weak and electromagnetic-weak currents and can be used for obtaining the corresponding reduced 
meson currents entering to the bosonized nonleptonic weak Lagrangians. 

In such approximation the problem of double counting does not arise. The effect of TrAi-mixing, 
being most important for the description of radiative weak decays, is taken into account by the 
corresponding TrAi-diagonalization factor. This factor appears explicitly in the bosonized strong 
Lagrangians and weak and electromagnetic-weak currents after reduction of the vector and axial- 
vector fields. With such modification of the strong Lagrangian and currents it is possible to reproduce 
within the nonlinear parameterization of chiral symmetry most of the results of the linear Lagrangian 
approach concerning ttAi -mixing effects and meson resonance exchange contributions. 

In Section 1 we discuss the basic formalism and display all definitions and constants related to the 
bosonization of quarks in NJL model. For convenience, all cumbersome expressions resulting from the 
heat-kernel computation of the quark determinant are given in the Appendix A. The total expressions 
for the bosonized effective Lagrangians including vector and axial-vector fields are presented in the 
Appendix B up to O(p^) terms. In Section 3 we consider the static equations of motion for chiral 
rotated collective meson fields in unitary gauge. Applying these equations of motion we eliminate the 
heavy meson resonances from the modulus of the quark determinant and obtain in such a way the 
effective pseudoscalar strong Lagrangian with reduced vector and axial-vector degrees of freedom. 
The reduced pseudoscalar {V — A) and {S — P) currents corresponding to the respective quark 
currents and quark densities are obtained in Section 4. In Section 5 we discuss the results of some 
numerical estimations and phenomenological analysis of the structure constants for the reduced strong 
Lagrangian and currents. 

2. Bosonization of the NJL model 

The starting point of our consideration is the NJL Lagrangian of the effective four-quark interaction 
which has the form fl^ : 

C-N.JL = liid - mo)q + dnt (1) 

with 

Here Gi and G2 are some universal coupling constants; mp — diag(rnQ,mQ, ...^rriQ) is the current 
quark mass matrix (summation over repeated indices is assumed) , and A° are the generators of the 
SU{n) fiavour group normalized according to tr A"A'' = 2Sab- Using a standard quark bosonization 
approach based on path integral techniques one can derive an effective meson action from the NJL 
Lagrangian (|l|). First one has to introduce collective fields for the scalar (S), pseudoscalar (P), vector 
{V) and axial-vector {A) colorless mesons associated to the following quark bilinears: 
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After substituting these expressions into Cnjl the interaction part of the Lagrangian is of Yukawa 
form. The part of Cnjl which is bihnear in the quark fields can be rewritten as 



C = qiDq 



with D being the Dirac operator: 

= [i(d + Ar) - ($ + mo)]PR + [i{d + Al) - ($+ + mo)]PL. 



(2) 



Here ^ = S + iP,V = V^7^, A = A^'y^'] Pl,r = ^(1 ± 75) are chiral projectors; Al,r = ± A are 
right and left combinations of fields, and 



s = s^-, 



\a \a \a 

P = P<^ y A^ = ~iAl — 



are the matrix- valued collective fields. 

After integration over quark fields the generating functional, corresponding to the effective action 
of the NJL model for collective meson fields, can be presented in the following form: 



Z = J V^V^+VVVA exp[i5($,$+,y, A)] , 



where 



1 



1 



— tr[(cl> - mo)+(<I> - mo)] - -^tr{V' + A') 



4Gi 



4G2 



(3) 



i Tr' [log(iD)] (4) 



is the effective action for scalar, pseudoscalar, vector and axial-vector mesons. The first term in (Q), 
quadratic in meson fields, arises from the linearization of the four-quark interaction. The second term 
is the quark determinant describing the interaction of mesons. The trace Tr' is to be understood as 
a space-time integration and a "normal" trace over Dirac, colour and fiavour indices: 



Tr' 



ct^x Tr , Tr 



tr. 



tr r ■ tr 



The quark determinant can be evaluated either by expansion in quark loops pO|-]T2{ or by the heat- 
kernel technique with proper-time regularization |^ . Then, the real part of log ( det iD) contributes 
to the non-anomalous part of the effective Lagrangian while the imaginary part of it gives the 
anomalous effective Lagrangian of Wess and Zumino which is related to chiral anomalies |Q . 

The modulus of the quark determinant is presented in the heat-kernel method as the expansion 
over the so-called Seeley-deWitt coefficients hk'- 



log I det iD I 



1 r{k-2^l/I^ 

,,2k 



2 (47r)- 



where 



T{a, x) 



(5) 
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is incomplete gamma function; fj, plays the role of some empirical mass scale parameter which will fix 
the regularization in the region of low momenta, and A is the intrinsic regularization cutoff parameter. 
There are some technical reasons for using the heat kernel method instead of the "straight" method 
of calculating quark loops. The main advantage of this method is that its recursive algorithms 
can be adopted on Computer Algebra Systems such as FORM or REDUCE quite effectively. In 
particular higher order derivative contributions can conveniently be calculated, too. The formulae 
for the Seeley-deWitt coefficients /i^ up to fc = 6 obtained in |29j and the full expressions for p'*- and 
p^-contributions to the bosonizcd meson Lagrangian are presented in the Appendices A,B. 

We will consider a nonlinear parameterization of chiral symmetry corresponding to the following 
representation of $: 

where T,{x) is the matrix of scalar fields belonging to the diagonal flavour group while matrix Q{x) 
represents the pseudoscalar degrees of freedom tf living in the coset space U{n)L x U{n)ii/Uv{n), 
which can be parameterized by the unitary matrix 

n{x) = exp i^--^ip{x)j , ip{x) = , 

with _Fo being the bare tt decay constant. Under chiral rotations 

the fields $ and ^^/^ are transformed as 
and 

For the unitary gauge — ^r — Vl the rotated Dirac operator (|^) gets the form 

{D->ii)^ {PLn + PrVl+)S{PlVl + PRn+) = ^(9 + ^ + 175) -e. (7) 

It is worth noting that under local UL{n) x Unin) transformations the modulus of the quark deter- 
minant is invariant, while the quadratic terms of V^, and the chiral anomaly do not respect this 
invariance. 

Note that there arises a quark condensate <qq>^ owing to the nonvanishing vacuum expectation 
value of the scalar meson field So realizing spontaneous breakdown of chiral symmetry. In fact, 
assuming approximate fiavor symmetry of the condensate and using the equation of motion 5*0 ~ 
~2Gi^J^q one gets 

<S'o>= -2Gi\ - <qq>= i2Gi \ - Tr {iD($ = fi,V = A = 0)}"^ . 
V n y 11 



4 



Thus, our mass scale fj, = <So> is determined by Gi and <qq> or (using the explicit expression 
for the condensate by a loop integral) by d, the cutoff A and the current mass mo. 

Taking into account the equations of motion for nonrotated scalar and pseudoscalar meson fields 
in nonlinear parameterization one can derive from and eqs.( ^ , ^^ of the Appendix B the following 
general expression of the effective meson Lagrangian including p^- and p^~interactions: 

47/''^ = -^tr{L^L^^) + ^tT{MU + U+M) 

Li - ^L2^ ( ti-L^Lf^y + L2 tr L,]^ + 3(L^L^)'^ + L3 tr {{L^L^')^) 

+Uir {D^Ul3^U+) tr M{U + U+) + L^iv D^U~D^U+ {MU + U+ M) 
+Lg (tr {MU + U+M)y + ^7 ( tr {MU - U+ M)^ + Lg tr {MUMU + U+MU+M)) 
+L9 tr (fI+^ D^'UD" U+ + fI^-W"U+ D-'U^ - tr (u+ F^+^UF^-'^'"'^ 
+Fltr((FW)%(F(;))^)+i^2trM^ (8) 

where the dimensionless structure constants Li{i — 1, 10) and Hi 2 were introduced by Gasser and 
Leutwyler in ref. ||2^. Here we have introduced the notations 

U^n^ ; L,^D,UU+ ; F,^ ^ , 

with y = r(0,AiVA') ; M = diagixl,xl-,xl) , X^^ ^ rnl^il [GiF^) ^ -2w'^<qq>F^^ ; 
< qq > is the quark condensate; 

are field-strength tensors, and 

D^* = d^*+{A^^*~*Al+^), D^* = d^*+{Al+K~*Al-')) (9) 

are the covariant derivatives; A^/^^ — zL A^. Moreover, the coefficients Li and Hi 2 are given by 
Li — ^L2 = L4 = Lq = and 

Nr 1 



L2 = 


Nc 1 
167r2 12 ' ' 


L5 = 


Idtt^ 


Ls = 


Nc 2 r 

Idtt^ 


Hi = 


Ne If 1\ 

16^2 6 I ^ 2 / ' 



167r2 6 ' 



167r2 6 12 

A^. 1 , _ AT, 1 
^10 — ^ 



167r2 3 ' 167r2 6 ' 



where x = — /^Fq2/(2 <qq>). 
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The effective (nonreduced) Lagrangian for the pseudoscalar sector, taking into account also the 
emission of the "structural" photons y^^, can be obtained from (||) when VJj = = and tensor 
-F/iJ'' is replaced by ie{dfj,A^ — d^Ai_i,). In the following section we will discuss the reduced nonlinear 
Lagrangian for pseudoscalar fields, which arises from generating functional (|^) after integrating out 
the vector and axial-vector degrees of freedom in the modulus of quark determinant. 



3. Strong Lagrangians with reduced vector and axial- vector 
fields 

To perform the integration over vector and axial-vector fields we will use the fact that the modulus of 
quark determinant is invariant under chiral rotations. Then, the pseudoscalar fields can be eliminated 
from the modulus of quark determinant in the effective action (jj) by using the rotated Dirac operator 
(|^) for unitary gauge. After such transformation the pseudoscalar degrees of freedom still remain in 
the mass term of eq.(|]), quadratic in meson fields, which are not invariant under chiral rotations. 
Since the masses of the vector and axial-vector mesons are large compared to the pion mass it is 
possible to integrate out the rotated fields Vp, and (^) in the effective meson action using the 
equations of motion which arise from the mass terms of the effective action (jj) in the static limit 
||30|] . In such an approximation the kinetic terms (i^^iJ'')^ for the rotated fields andA^ as well as 
higher order derivative nonanomalous and Wess-Zumino terms are treated as a perturbation. 

In terms of the rotated fields V^, @ the quadratic part of the effective action (||) leads to the 
Lagrangian 

£o = ^ tr {MU + h.c.) - i^-fj tr [{V,. - v^f + {A^ - a^f] , (11) 

where {rny/gy)^ — 1/(402), with rriy and gy being the bare mass and coupling constant of the 
vector gauge field, and 



The modulus of quark determinant contributes to divergent and finite parts of the effective meson 
Lagrangian. In terms of the rotated fields, taking into account that for unitary gauge <& ^ S, the 
divergent part of the quark determinant (eg. (pi]) from the Appendix B) gives 

- ||tr{i^,I]i5^I]+ + i[(FW)%(F(;))^]} 

= lltr!^[V,,mor-{A,,mor-M^rno + fi)Al + ^[{F(t^^^ (12) 

where the approximation S = /i + mo was used. 

The p'^-terms of the finite part of the effective meson Lagrangians (eq.(p^ from the Appendix B) 
are of the form 



Ap I _ tr ^ - 



Nr ( 1 



1 



6 



3 V 



(-)\2 



tr 



2M 



,({A.,{A^mo}}-[K,[^^mo]]) 
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(13) 



where toq = mo//i, and we used the approximation A1 = = — /i^ « for matrices 

and in unitary gauge. 

The kinetic terms (F^lf'"^'), arising from the sum of Lagrangians (^2|) and (|l3|), lead to the standard 
form after rescahng the rotated nonphysical vector and axial-vector fields V^, A^: 

Here 



r _ 9V y{ph) J _ 9V J(ph) 

(1+^)1/2 V ' (^_~)1/2^M 



(14) 



487r2 V NeH^ 



1 



1-1/2 



487r2 
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and Vji^'*'' , vl^^'*'' are the physical fields of vector and axial-vector mesons with masses 



\2 



ml = 



1 + 7 



(15) 



(16) 



where = 1 — {Fogy /rriy)'^ is the TrAi-mixing factor. 

Since in the following we also want to investigate the radiative processes with "structural" pho- 
ton emission in addition to inner bremsstrahlung ones, it is necessary to include electromagnetic 
interactions in the bosonization procedure. Obviously, one then simply has to use the replacements 

, or — ^ + ieoAp,Q , where Q is the matrix of electric quark charges 

and 

(1+7)1/2 



Aiph) _ 9v . (ph) _ 



9'v 



are the physical electromagnetic field and charge respectively. 

The static equations of motion arise from variation the mass terms of eq.(|ll|) in chiral limit over 
rotated fields V^, and lead to the relations 



V = w^'*'' A = Zl a^'^^ 



(17) 
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and 



?(±) 



(18) 



Here 



,(7) _ 



d\jj'^* — d^*+ie^A^[Q, *] = d^*+ie^P'^'>A''iI''^\Q, *] is the prolonged derivative describing the emission 
of the inner bremsstrahlung photon while the electromagnetic field strength tensor TjZ' = d^Av — 
diyAp, corresponds to the structural photon (ei^T^jj) = e^^'^^ tIjI>^^^)', and L^^^ — {d']jJ^U)U^ . Further, 
we will omit for simplicity the upper indices (7) corresponding to the inner bremsstrahlung photon 
and only tensors J-\jl) will be kept explicitly. We will also omit everywhere the upper indices {ph) 
assuming that all photons and electromagnetic charges in further formulae arc physical. 

Applying the equations of motion to the terms of the effective actions (|ll|,^2|), quadratic in 
vector and axial-vector fields, one reproduces the standard kinetic term for the pseudoscalar sector: 



f tr(L,i-). 



(19) 



In the same way the p^'-terms of the actions (12 ,1^) lead to the reduced Lagrangians for pseudoscalar 
mesons of the types 



^Z^tr 



([i^, L,]2) + (3L2 + is) tr ((i^X^)' 
2Z5tr (L^L'^CflA/eA) - 2ie.F^^)Z9tr {q^+L^^L-(,i 



2(ze)^Lio tr Ai Q^+L,^r Q^W^r ~ Q'^IK^R 



(20) 



n(p\red) ^ J^r. 4: 
32 



4 r 



1? 



-^e^Z^^i tr {q[(1 + Z\){C+L^^L^^R, mo} + 'iZ\i^+L^^Rmoi+L-'iR)\ } 

-(*e)2^izl tr[mo[Al(2iQC+L,(Rf^^+L,^R{i+L,^R,Q'}) 

-A,Au {2Qi+L^tR Qi^L'^^R - i+L^tR&'^^R, Q2})] I + o{rnl) , (21) 

Here C^'p '^'^''■^ represents the part, corresponding to the effective p^-Lagrangian in the Gasser- 
Leutwyler representation with the structure coefficients Li defined by the relations. 



IGtt^ 



lzui(.l-i)((^l-i)^i^-^.i 
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= 

Lg = 



167r2 



167r2 



-^^10 — ^Zj^ - 



1 1 + 7 



(22) 



where 



1+7 _ 



The Lagrangian Cm„'^'^'^'' describes the additional corrections arising from the expansion over the 
quark mass toq. 

To take into account also the p^~terms of the finite part of the effective meson Lagrangian (eq. ( ^ ) 
from Appendix B) we will restrict ourselves only by consideration of the terms which additionally 
contribute to TOQ-corrections of type of eq. (pl|) : 



-fin 



Nr 



327rV* 



■tr 



M 



M{{D^,YD^Y+Y - {Df,Y D^Y+y + D^,Y D^Y+ D-'Y d'"Y+ 



12^2 



{ D^Y+D^YM + M Z?^E+ 0"^ + D^Y+ M D^Y 



M{Flz^f + M{Fl+^f] + (other terms) , 



(23) 



All other omitted terms after reducing the vector and axial- vector fields in eq.(^3|) will contribute 
only to the interaction of six or more pseudoscalar mesons. 

The final expression for the reduced Lagrangian Cmg'^'^'^'^ can be presented then in the general 
form 



'-'mo 



(lefQ^ tr {mo [^^ (2(Q^jL.Cfl)' - ilL,iB{i\Uin, Q^} 



(24) 



with Qi being the structure coefficients: 



Qi = 



Nr. 1 



167r2 6 
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1 

167r2 3 

Nc 1 
167r2 6 

A^c 1 
167r2 12 



16Z1 



-(Zl-l)(17Zl-9) 



16 



3Z% + 16ZI 



7), Qs = 0, Q, 



-(Zl-l)(7Zl-9) 
Nc 3 



74 

167r2 8 



(25) 



4. Reduced currents 



The path-integral bosonization method can be applied to the weak and electromagnetic-weak currents 
by using a generating functional for Green functions of quark currents introduced in p9| and |20 



After transition to collective fields in such a generating functional the latter is determined by the 
analog of formula (||) where now iT) is replaced by 



iBii]) = [iid + Ab.- ivn) - ($ + mo - Vb)]Pr^ + [iid + Al - ifjn) - ($+ + mo - vl)]Pl 



(26) 



Here rjL,R = Vl fi/i7^^ ^^'^ Vl.r — vt _r V ^^'^ external sources coupling to the quark currents 
qPl,r'J'^ a-iid quark densities qPh^R^Q respectively. The quark densities define the contributions 
of the penguin-type four-quark operators of the effective nonleptonic weak Lagrangian 111 to the 
matrix elements of relevant kaon decays. The bosonized {V =F A) and {S ^ P) meson currents, 
corresponding to the quark currents qPL.RliJ.\-<l a-nd quark densities qPlm^-Q: can be obtained by 
variating the quark determinant with redefined Dirac operator ( p6| ) over the external sources coupling 
with these quark bilinears . 

For further discussions it is convenient to present the bosonized weak and electromagnetic-weak 
{V — A)-current for pseudoscalar sector, generated by the nonreduced Lagrangian (|[) and including 
the electromagnetic-weak structural photon emission, in the form: 



J 



nredia 



^^tr (a% 
1 



-i tr ■ 



-RiU+{iMU + C/+M), L^}C/ + i?2L.L^i" 



+eJ^jZ^R,ty(^y^[UQU+,L'' 



(27) 

Here, the first term is the kinetic current and all other terms originate from the p*-part of La- 
grangian (H); Ri are the structure coefficients, associated with the corresponding parameters Ri of 
the representation (BOf) for the reduced {V — A)-currents: 



i?i — — L5 



i?2 — 2L2 , ^3 — 2L2 + L3 



1 

Ri — ^2^9 : 



Rf, — La 



^10 



The bosonized {S — P) current for pseudoscalar sector, generated by the Lagrangian 
including the structural photons, has the form: 



(28) 
]) and 



(nred)a 



= tr {X-d^U) + ^^lRtv {X^U) 
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1 

2 

ie 



i tr {a° [L2d^{L,L''L'') + {21^ + is)^^ {L.L'^L^) 



-Lg tr 



(29) 



Here, the first and second terms are generated at p^-level by the kinetic and mass terms of Lagrangian 
(^), respectively, while all other terms originate from its p^-part. 

Combining the method of the chiral bosonization of quark currents with the static equations 
of motions it is possible to obtain the bosonized meson currents for pseudoscalar sector with the 
reduced vector and axial-vector degrees of freedom. In this way one can reproduce the standard 
kinetic {V — A) current for pseudoscalar mesons 



J 



{kin)a 



which arises from the terms of effective actions (|11|,|^ , quadratic in vector and axial- vector rotated 
fields, after redefinition of the rotated fields 



V„. 



Vf, - i{r]Lf, + VB.ti) , Af,^ Af,+ i{T]Lf, - , 



and variation over ry/,^ with applying the static equations of motion. 

Applying the same procedure to the p^- and p^-terms ( [l3| , p3|) of the effective action we also 
obtain the bosonized weak and electromagnetic-weak (V — A) currents for pseudoscalar sector with 
the reduced vector and axial-vector degrees of freedom. It is convenient to present these reduced 
currents in the form: 



j: 



{p ,red)a 



j: 



{p ,red)a 



-ztr{A'^[i?2(eAi.iM^"Ci?.) 



-itr|A° 

hi?8 



(30) 



1^.. ^&{mQ,S,^Li,L^,L^^n} + RT{m^S,^LlL^,^R + ^^L 
+Rw^R(^L^L^(,Rmo^RL'' + L^^RTfio^RLf.L.y^^R 

+Riidu{il{L-^Rmoi+L^ - L^^Rmo^+L'')^R) + Ri2{rno, d,{^+[L^, L'']^r)}'\ } 



^Ru{moQ^+L''^R - i+L^iRQma) + Ri5{Qmo^+L-'^R - CA^'efl^oQ)] } 



(31) 
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with Ri being the structure coefficients 

Nc 1 o , 

^1 = -TT^oZUy-^), 



167r2 2 



2 ^ 167r2 24 (5y)^ 



167r2 6 "-^V^ {cfvY 
^ 167r2 96 ^6 167r2 24 

^" = 16^^12^^^^^"^)' 



Rl3 = ^^14 = ii?15--^^^l. (32) 

Thus, the reduction of the vector and axial- vector fields does not change the kinetic term of the 
bosonized {V — A) current while the structure of the p^'-part of {V — A) current is strongly modified 
(compare ( p7|) and (|3C|)). 

Using the bosonization procedure of ref. j2^ and the equations of motion (|l7| ) we obtain also the 
reduced {S — P) meson currents. After redefinition of scalar fields 

S ^ E - 277i^ , S+ ^ S+ - 27jl (33) 

and variation over r]L with applying the static equations of motion the divergent part of the effective 
action (|l^ ) leads to the scalar current 

^f "'^^'^" = tr { A'' [^Ll^n + V^o + \ (2(^+L,^^moCAi''e«) + {mo, CA^^CW)] } 

+ -f^/xi?Z^'tr(A"(l + mo)), (34) 



where R —<qq>l(^\xFl). 

The "p^- and p^-terms of the finite part of the effective action (^^3|) generate the scalar meson 
current which can be present in the form 



+G4(a^(e+L.?fl,) eAi^i^Cfl + eA^'i^Ci? a^(eAi.Ci?, 
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+G^{eRLlLUu) + G^{£+{L^L,f£,R)\ } 

+ {ieT^;Jf-Giz tr (A'^Q^^ + 0(mo) (35) 



where G,- are the structure coefficients: 



G, - 6G2-6G3~^Gn-Y|^Zl, ^4 - G5 - , 

Nr. 1 



- 256^ 6 

= 256^^3^^' ^^^ = " 256^^^^^^"^)' ^^^^ 256^3- ^^^^ 
It is important to mention that the p^-terms (|2^) of the finite part of the effective action contributes 
only to the mo-corrections ( ^ ) for {V — A) current while in the case of {S — P) currents due to the 
fact that redefinition (|33| ) leads to the replacement 

the same p^-terms give contributions also to the part of current ( |35| ) do not containing quark mass 
mo. 

It can be easily shown that at the p^-level the reduction of the vector and axial-vector fields 
does not change the physical results for matrix elements of the bosonized gluonic penguin operator. 
In fact, both for the reduced current (^5|) and for the first two terms of the nonreduced current 
( |29| ) the corresponding -contributions to the penguin operator matrix element, dominating by the 
interference of the kinetic and mass terms of the scalar current, can be presented effectively in the 
same form: ^ 

< >^ _ ^ {d^Ud^U+)23 > ■ 

On the other hand the structure of the p^-part of the pseudoscalar meson (5* — P) current proves to 
be strongly modified by the reduction of the vector and axial-vector fields (compare the expressions 
(H) and (H)). 
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5. Numerical estimations 



To discuss some physical consequences for pseudoscalar nonet of mesons we have to fix initially 
the numerical values of the various parameters entering in the reduced Lagrangian and currents. 
The parameters xf can be fixed by the spectrum of pseudoscalar mesons. Here we use the values 
xl = 0.0114GeV2, xj = 0.025GeV2, and Xs = 0.47GeV2. To fix other empirical constants of our 
model we will use the experimental parameters, listed in Table 1: the masses of p- and ^i-mesons, 
the coupling constants of p ^ tttt and Tr,K ^ jiv decays, the TTTr-scattering lengths af, the pion 
electromagnetic squared radii < and pion polarizability a^ri . We also include in our analysis 

the data on the 77 — > tt+tt^ cross section near to the threshold (see Fig.l). We will use the relations 
(|l|), (|ll),5y = 5^(l + 7)-'/' and 



SpTTTT — gV 



1 



2Fi 



f2 
^0 



487r2 24^2 



Then, using the connection = mox^/m^, where mo 
F-n K+ can be presented in the form: 



(m° + m") /2, the splitting of the constants 



= Z. 



l+(mO+mylzl(^l 



™0 



:Zl(l 



izi 



1 - 



5Zl 



Here 



1 



{ml + to2,. 



|4Zl 1 



Wn^Fg 



1/2 



is the factor arising from the renormalization of the pseudoscalar meson fields 



1 



167r2i^2 



mo 



which leads the bilinear kinetic part, including terms up to order m^, to the standard form of the 
kinetic Lagrangian: 



rbii _ t tr 

'-kin ~ 2 



1 - 1 - 



16tt^F^ 



Too 



■ tr (9^$(9'^$) 



The TTTT-scattering lengths are defined by the structure coefhcients L2 and L^. For TTTr-scattering 
lengths aj (indices / and I refer here to the isotopic spin and orbital momentum, respectively) in 
one-loop approximation we obtained p2] 



-aoi9 - 56) + -al 



-ao2S 



" 9„ 



A 



5A + 3B + 2D + 3C - 6(S^ + 46 
-I- £» - 3((J2 + 6 + 3) 



3) 



14 



2«o+2«03 



B + 6S + a-3 + ^{d^ -b-3) 



4 



2«o 



2«o 



3(5 - 2a + 6 (52+46 + 3 

^ ^ 

65 - a + 3 2 
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9 



^((5^ + 6 + 3) 



Here ao = ^ (m^/(27rFo)) ; S — |(1 — /?), with /3 being the parameter of chiral symmetry breaking 
which takes here the value /3 = 1/2 ; a = 21(1 — 6) and b = (11(5^ — 15(5 + 3). The parameters 

A = A^+ A'°''P , B = B^ + B^°°P , C ^C^ + C'°°p , D = + D^°°p 

include in themselves the Born contributions 

= -1447r2(Z2 - Z3) , B^ = -5767r%, C^^ = 5767r2(Z2 + Z3) , D'^ = bl&'K'^L2 

and the pion-loop contributions calculated, using the superpropagator method Jssf , in ref . | |3^ : 



Table 1: Physical input parameters used for the fixing of the empirical constants of the model 



Input parameters 


Experiment 


Theory 


■niAi 
Fk 


770Mey 
1260Mey 

6.3 
93Mey 
119Mey 


772A/ey 
lieOAfey 

6.8 
93.9Mey 
\\%.hMeV 


al ■ m.^ 
al ■ rrij, 
ai ■ 

a'^-ml 
al-m% 


0.23 ±0.05 |35|1 
-0.05 + 0.03 H 
0.036 + 0.010 H 
(17 ±3) • 10-'' |6) 
(1.3 + 3) • lO-'^ Hi 


0.20 
-0.04 

0.038 
17-10-4 
2 • 10-'' 


< > + 

em TT^ 
"7r± 


(0.439 ±0.030)/m2 |37| 
(6.8 ± 1.4) • lO-Vm^ m 


0.53/m2 
8.0 • 10-''/m3 



The pion electromagnetic squared radius is defined as the coefficient of the (/■^-expansion of the 
electromagnetic form factor /^'"(g^): 

< i^{v2)\yr\<vi) /r('z')(pi -p2)m, 

/r('z') = l + J<'-em>. 9' + - 
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Being restricted only by pion loops, one gets in the SP-regularization the corresponding contribution 
to the electromagnetic squared radius 



^ ' em ^7T+ 



(47r^^o)2 



3C + ln 



V 27rPo ' 



0.062/m2 



where C = 0.577 is the Euler constant. Because the main contribution to this value arises from 
the logarithmic term, the kaon loop contribution, containing the small logarithm In ^mx/(27r_Fo)^ , 
can be neglected. At the Born level, the contribution to the pion electromagnetic squared radius 
originates from the Lg-term of the reduced Lagrangian (|2^) : 



(Born] 



< '"em >7r+ 



^0 



The pion polarizability can be determined through the Compton-scattering amplitude: 



< 7ri(pi)7r2(p2)|S'|7Ai((7i)7A2(92) >= Ti(piP2 | qiq2) + T2{piP2\qiq2) , 



P1P2 _ P1P2 

Piqi P2qi 



) = 0; 



T2=e- 



,£x^{{qiq2)g,,u ~ qii^q2,,)(3{qiq2), 



where P{qiq2) is the so-called dynamical polarizability function. Defining the polarizability of a 
meson as the coefficient of the effective interaction with the external electromagnetic field 



one obtains 



/37r(gig2) 

87rm7r 



(?ig2)=o 

The pion-loops give the finite contributions without C/y-divergences [6,7,11]: 



o(loop) 



where = {qiq2)/{2ml) , /(C) = J\0 - 1 , and 

arctan(C-i - 1)^1/2 

^(0 = 



[loop] 



Mn ^+V^EEL 
2 -i+yi^ 



47r2F2V^ 3s-. J 



0<C< 1; 
1; 
C<0. 



The meson-loop contributions to the pion polarizabilities are 

^2 



a^'T"^ = 0, 



(loop) 



3847r3F2m2 



-5.43 • IQ-^fnf 
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At the Born level, the ig- and Lio-terms of the reduced Lagrangian ( pO| ) give: 

o 2 _ _ 

o(Born) _ 2'_/r I r \ aiBorn) _ „ 

^0 

In our analysis the constants Fq, toq, ^J^ and my are treated as the independent empirical param- 
eters and their values are fixed as 

Fo = 91.7MeV, mo = 2.2MeV, /i = 186MeV, TO^= 840MeV. (37) 

The corresponding calculated values of the input parameters are also presented in Table 1. The 
results for the 77 tt+tt^ cross sections are shown in Fig.l. All other constants can be calculated 
using the values (|3^): 

glr = 5A, 7 = 0.185, Zl = 0.653, <gg>= -(330MeV)^ TO2 = 53MeV. 

The values for the current quark masses seem to be by a factor of 2 -f- 3 too small and the quark 
condensate, respectively, by the same factor too large as compared with the corresponding value from 
the usual phcnomcnological analysis, based on nonreduced Lagrangian and currents. A similar shift 
of the current masses and the condensate was already observed, for example, in ref.|Q after taking 
into account the vector-scalar and axial-vector-pseudoscalar mixing in the analysis of the collective 
mesons mass spectrum within the extended NJL model. 

Using the values of the parameters Z'j^,j and (5^)^ which were fixed above, one can compare 
numerically the structural parameters Li (^2|) of the reduced effective Lagrangian (|20| ) with the 
corresponding parameters Li of the nonreduced Lagrangian (j^): 

Z2 = 1.20X2 = 1.90 • 10-3, Z3 = 1.71X3 = -5.41- 10-3, Z5 = 1.99- 10-3, 

Zg = 1.35X9 = 8.53- 10-3, Zio = 1.36X10 = -4.33- 10-3, (38) 



After substituting the values of ,7 and (gy)^ into eqs.(p2|) one can also compare numerically the 
structure parameters Ri and i?^: 

Ri = -0.285 • 10-3 , R2 = 0.76i?2 = 2.42 • 10-3 , ^3 = -0.992 • 10-3 {R^ = 0) , 

^4 = 2^3 = 0.62i?4 = -1.98 • 10-3 , ^5 = 0.39i?5 = 1.23 • 10-3 . (39) 



The electromagnetic-weak part of the non reduced current ( |27D corresponding to the structural 
constant R5 (respectively, the i?5-term of the reduced current (^) describes the axial- vector form 
factor Fa of the radiative decay tt — > Ivj. The form factors of this decay are defined by the parame- 
terization of the amplitude 

T^{K,TT~^li^-f) =\/2e[Fv£^,a0rg"e'' + zFA(£^(fc(7) -(7^(fc£))], 

where k is the 4-momentum of the decaying meson, q and £ are the 4-momentum and polarization 
4-vector of the photon, and the vector form factor Fy is determined by the anomalous Wess-Zumino 
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electromagnetic-weak current, originating from the anomalous part of the effective meson action, 
which is related to the phase of the quark determinant. The ratio of the axial-vector and vector form 
factors is determined by the relation 

Fa 9 
tv 

The theoretical value of the ratio Fa/ Fy — 327r^(Lg + -Lio) = 1 arising from non reduced current ( ]27| ) 
with structure constants Lg^io ( p^ is in disagreement with the experimental results on this ratio: 

[exp] 



^^yexw ^ 1 0.25 ±0.12 
~ 1 0.41 ±0.23 



At the same time the i?5 gives the value 



Fa 9 / 127r^ 1 + 7 \ 

^ =-Z5 1--— =0.39 



Fv {g'vY 

in agreement with the experimental data. 

This result reproduces in nonlinear realization of chiral symmetry the well known result of ref.| 13 
where the role of TrAi-mixing in tt — > /z/7 decay was considered in linear parameterization. Thus, 
after reducing the vector and axial-vector degrees of freedom it proves to be possible to remove the 
inselfconsistency in the description of the ratio Fa/ Fy and pion polarizability which arises seemingly 
in the pseudoscalar sector of the non reduced effective Lagrangian (||) with the Lg^io-terms [| (see 
the detailed discussion of this inselfconsistency, for example, in ref. ^). The same problem 

was also considered in ref.[Q, where the values of the structure constants combination (Lg + Lio) 
and pion polarizability a^i determined from the fit of 77 —t tt+tt^ cross section data were discussed. 
Fig.l shows that within the experimental errors the MARK-II data |^ are consistent with the 
experimental result for pion polarizability obtained from radiative tt scattering in nuclear Coulomb 
fields |38j. We have taken into account one-loop corrections, while this was not done in ref.[Q. The 
description of the 77 tt+tt^ cross section data above m^rTr — 500MeV can be improved if one takes 
into account the unitary corrections in a more complete way pffl. 



Conclusion 

In this paper we considered the modification of the bosonized Lagrangian and of the currents for 
the pseudoscalar sector obtained after integrating out the vector and axial-vector collective fields in 
the generating functional of the NJL model. It has been shown, that the reduction of the meson 
resonances does not affect the kinetic terms of the strong Lagrangian and the bosonized {V — A) 
current as well as the (5' — P) current, generated by the divergent part of quark determinant. On 
the other hand, the reduction of the vector and axial-vector fields leads to an essential modification 
of those part of the pseudoscalar strong Lagrangian and of the currents, which originate from O(p^) 
terms of the quark determinant. The reduced Lagrangians and currents allow us to take into account 

^This point was already noted shortly in our resent paper ]44| . 



18 



in a simple way all effects arising from resonance exchange contributions and Tr^i-mixing when cal- 
culating the amplitudes of various processes with pseudoscalar mesons in the initial and final states. 
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DFG for support of investigations (Project Eb 139/1-1). Another of the authors (A.V.Lanyov) is 
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Appendix A. Heat-kernel computation of quark determinant 

The logarithm of the modulus of the quark determinant is defined in "proper-time" regularization as 

1 



^ ^ 1 /■ 1 ^ ^ 

log|detiD| = - Tr' log(D+D) = - - / dr - Tr' exp( - D+Dr) 
2 ^hih? T 



(40) 



with A being the intrinsic regularization parameter. The main idea of the heat-kernel method is to 
evaluate 

< X I exp(-D+DT) I y > 

around its nonperturbated part 



< X I exp(-(n + ^2)r) I y >= 



(47rT)2 

in powers of proper-time r with the so-called Seeley-deWitt coefficients hk{x, y) 



< X I exp(-D+Dr) | y > = 



(47rT)2 



-M'T + (x-y)V(4 



After integration over t in (|4^) one gets the expression (^) for log | det zD | . Using the definition of 
the gamma function r(a, x), one can separate the divergent and finite parts of the quark determinant 



1 



- log(det D+D) = Bpoi + Bjog + B^^, 



where 



B„„^ — 



2 

1 e-^ 



Pol 2 (47r)2 



—T Tr' ho + -(— Tr' ho - fi'^ Tr' hi) 
Ix'^ X 1 



has a pole at a; = ^ / — 



^iog = -^(i^r(o,^) 



— /i* Tr' ho — fi^ Tr' hi + Tr' /i2 



is logarithmic divergent, and the finite part has the form 



%n - 



1 1 



2 (47r)- 



^ /-2'=r(A: - 2, a;) Tr' /ifc. 



fc=3 



The very lengthy calculations of the Seeley-deWitt coefficients hk can be performed only by 
computer support. The calculation of the heat-coefficients is a recursive process which can be done 
by Computer Algebra Systems such as FORM and REDUCE very conveniently. In ref.||2^ we have 
calculated the full coefficients up to the order fc = 6. After voluminous computations one gets the 
complex expressions for heat-coefficients hi, . . . h^: 



ho{x) 
hi{x) 



1, 
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+ 1^^^'^^""'^^''' " 12^ (29r'^„r^" + 27r^"r^) (x;^, + i^,^^) 

(OO A 1 27 1 \ 

25200 + 1575 5040 600 



1 Q 47 17 

12600 « 16800^ ' ' ' ' ' 25200^ 



1575^ ' 25200 12600^ ' 1575^ 

6300 ^ 5600^ '^^ 5040 ^ ^ 



Here, 



Sjj, = [dfj,,a], Sjj,u = [d/j,, S,^], Sxfj,,^ = [dxjS^i,] , etc. 
are commutators of the operators and a which are defined by the relations 

d^ = 5^ + r^, r^ = v^, + A^j'', a(x) = ^vi^ + ^^+i^ + i[7^71V-M'• 

We used the following notations: 

= Pfl$ + Pl^+ = s + , r^, = [d^, d,] = d^T, - d,T^ + [T^, r,] = F^^ + i^F^t 
with -F^"^ as field strength tensors, 

F^, = d^A, - d,A^ + [V^, A,] + [A^, K], 
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and 

V^H = d^H + [V^, H] - 7^{A^, H} 

as the covariant derivative. 

In the same way the next orders of heat expansion coefficients hi can be obtained using the 
developed computational technique based on the usage of computer algebra |2g[| . For simplicity we 
present below expressions only for minimal parts of heat-coefficients, i.e., only for the parts which do 
not vanish in the pseudoscalar region of the theory when = Af^ = 0: 

11 IQ 2 

+ 1008^'^'^ + 2800*^^"^* + 225*'^*"* 



25200 



1 



20160 



25200 ^ '^ ' *^ ' ^ ' 2016 

1 „ „ 1 



(375^'' + 235"''') + z—:aSf, {S^'%°'c. + 5 



9450 ''^ ^ 9072 ^ 

'"^^ (isSo^^^"' + sSo^"^'^'^ + 10^^'^'^'') + 252' 



-aS^fj^i, {52Sa°"^ + 535""^^) + aS^^ i, I 5"^^ TToTTi'^Ma" 



352800 /-"^ " " V3360 ^ 11340 



226800 ' ' 317520 

[53(5^5'')2 + (5^5.)^] + ^^S,S. (1575''„"'' + 2985'"'"„) 



70560 ^ " " ■ " ' 720 

5 1 

14112'^'^'^'"^""^'^ + K)5840'^'^ [5''"' (37S'"Q.,y + 7G>S'yc[") + 355'^'' S'",yc[] 

1 



21168 ' " 2880 



^^^^"'^ ( 141120'^''"" + 9^^""'' j + Soiio^^ ^^^^""^ + ^^^"""-^ ^'"^ 
+^5^ [2 (5""" + 5""" + 5""") + (5,„" + 25"<,,) 5""] + ^^^^.^.S""" 

^ (701(5/)3 + 5835'"' 5^„5^) - -^5"„(5'"')2 



83 1 „ 1 



635040 ^ ^ / ■ 316386 
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2835 952560 
'226800^ ''""^^ 12700800^ ^ 66150 '^"^ 

1 IS I 



52920 604800 

To obtain these expressions for the heat coefficients, we have made extensive use of the cychc prop- 
erties of the trace and of the Jacobi identity. 



Appendix B. Bosonized effective Lagrangians 



The effective meson Lagrangians in terms of collective fields can be obtained from the quark determi- 
nant after calculating in tr'ft,i(a;) the trace over Dirac indices. The "divergent" part of the effective 
meson Lagrangian is defined by the coefficients /iq, hi and /12 of the expansion (||): 



N, 



div 



327r2 
+2 



t'^^rio,^ 



M 



(41) 



where D'* and are covariant derivatives defined by eq. 
Applying the properties of the covariant derivatives 



and M = - 



D^{Oi02) = iD^Oi)02 + OiiD^02) = {D^Oi)02 + OiiD^02), 
D^{Oi02) = (D^0i)02 + OiiD'^02) = (d'^0i)02 + Oi(D^02); 



with 



+[4-),*], D * = a^*+[4+),*], 



and assuming the approximation V{k,ii^/K'^) k, V{k) (valid for fc > 1, and /i^/A^ <c 1) one can 
present the p''-terms of the finite part of the effective meson Lagrangian in the form 



fin 



N, 



3277^4 



■tr 



3 V 



(42) 



where 7W = <&+$- 
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In an analogous way, the p^-terms of the finite part of the effective meson Lagrangian can be 
presented as 



■'fin — 



327r2/i< 



tr 



1 r 

45 



60 



1 r 

'90 



1 

12 



180/i 

-3(15^$+£>'^$ £>^*+I?a$ + :d^$+D^$ iD^^+i)"* iD^^+D"^*) 



6^ 
1 , 

60* 



{D^M)' + {D^My 

^ V 5 



60 

180'^ 360^ 



- D''^D^<i>+) + Fj^+'>F^+'>''{D'"^+D''<^ - D^'i'+D"^) 
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Figure Caption 



Fig.l. MARK-II cross section data for 77 tt+tt^ for CMS production angles \cos9\ < 0.6. 
The experimental points in the region mTTTr <0.5 GeV were only included in the analysis. The dotted 
line shows the QED Born contribution; the dashed and dash-dotted lines show the results of the 
successive inclusion of ^^-contributions and one-loop corrections. Both lines are calculated with 
{Lg + Lio) = 4.2 • 10^^, corresponding to the fit of the total cross section data together with the 
parameters of Table 1. The solid line corresponds to the direct fit of the experimental points for 
fnjrTT <0.5 GeV without including the experimental parameters of Table 1. 
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